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A particle-hole symmetric Anderson impurity model with a metallic host of narrow bandwidth is 
studied within the framework of the local moment approach. The resultant single-particle spectra are 
compared to unrestricted Hartree-Fock, second order perturbation theory about the noninteracting 
limit, and Lanczos spectra by Hofstetter and Kehrein. Rather accurate analytical results explain 
the spectral evolution over almost the entire range of interactions. These encompass, in particular, a 
rationale for the four-peak structure observed in the low-energy sector of the Lanczos spectra in the 
moderate-coupling regime. In weak coupling, the spectral evolution is governed by orbital effects, 
while in the strong coupling Kondo limit, the model is shown to connect smoothly to the generic 
Anderson impurity with a flat and infinitely wide hybridization band. 
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I. INTRODUCTION 

Five decades of intense experimental and theoreti- 
cal research have boosted the Anderson impurity model 
(AIM)^ far beyond the scope of the "localized mag- 
netic states in metals" that it was initially designed 
for. In the first three decades, its concept of a sin- 
gle level with on-site Coulomb repulsion coupled to a 
host without electronic interactions was mainly used to 
describe magnetic transition metal impurities dissolved 
in otherwise nonmagnetic bulk metals. The last two 
decades' extraordinary progress in nanotechnology, how- 
ever, brought a myriad of new and rather surprising 
implementations of what in the meantime had become 
one of the th eoris ts' favourite toys, ranging from tunable 
quantum dot^^IH over carbon nanotube^^ and adsorbed 
organic molecule^ to single-electrorP or single-molecule 
transistors.'^ 

The vast majority of theoretical work, comprehensively 
reviewed in Ref. 9, focuses on AIMs with metallic hosts 
of large bandwidth. In the limit of infinite bandwidth, 
the exact static and thermodynamic properties can be 
deduced from the Bethe ansatz solution.'iSHIll A power- 
ful and versatile alternative is provided by the numerical 
renormalization group (NRG)p211Il recent extensions of 
which have also been able to address the dynamics of the 
model to excellent accuracy.fiSMSi] 

Yet another field of application is the Mott metal-to- 
insulator transition in high spatial dimensions PsHllliiej-e, 
dynamical mean-field theory (DMFTjPS! reduces the at 
first sight unrelated problem of interacting electrons on 
a high-dimensional lattice to an effective AIM immersed 
in a bath of identical sites whose properties have to be 
determined self-consistently. It was within this context 
that the possibility of an AIM with a metallic host of nar- 
row bandwidth was first evokecF^ since it naturally arises, 
in the vicinity of the Mott-Hubbard transition,'^ within 



the now widely accepted scenario of a metallic state sur- 
rounded by a preformed gap.f^^l At the time, W. Hofstet- 
ter and S. Kehrein argued, on grounds of their Lanczos- 
determined single-particle spectra for a corresponding 
AIM,'2SI that this incipient gap might be populated by 
localized states. Although these states have not actually 
been observed inside the preformed gap of the infinite- 
dimensional Hubbard model, recent high-resolution dy- 
namic density-matrix renormalization group (DDMRG) 
datcp21l22l Jo indeed show very narrow features on the 
inner band edges of the Hubbard satellites in the appro- 
priate regime of interactions. 

Independent from the relevance for the Mott-Hubbard 
transition, several questions about these sharp features 
may arise: (i) what are the underlying physical pro- 
cesses? (ii) do these processes depend on correlations 
within the DMFT bath? (iii) why do these features dis- 
appear for both, small and large values of the Coulomb 
repulsion? (iv) are they related - and if yes, in which 
manner - to the series of peaks observed in the low-energy 
sector of Lanczos-determined spectra of an AIM with a 
correlationless narrow-band host'PSl and (v) why, in the 
latter case, are the peaks organized in a four-set struc- 
ture? 

Some of the above questions have already been ad- 
dressed in a previous articlej^SSl in which the AIM was 
mainly studied perturbatively about the noninteracting 
limit; others, in particular those concerning the rich low- 
energy structure at moderate interaction strengths or the 
strong-coupling Kondo regime, lie out of reach for per- 
turbative approaches and remain unanswered thus far. It 
is these questions, among others, that shall be addressed 
in the present article, primarily within the framework of 
the so-called local moment approach (LMA). This non- 
perturbative many-body Green funct ion form alism, de- 
veloped by D. Logan and co-workers P^ * ^^ * ^ ^ introduces 
the concept of local moments from the outset. At pure 
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mean-field level, this would lead to a doubly degener- 
ate ground-state, as appropriate for an insulator, but in 
manifest contradiction to the Kondo singlet observed for 
impurities hosted in metals. The LMA aims to tran- 
scend this deficiency by accounting for dynamical tunnel- 
ing processes between the two mean-field ground-states, 
at a rate which has to be determined in consistency with 
Fermi-liquid behaviour on the lowest energy scale. The 
resulting formalism is equally well adapted to impuri- 
ties in metallic^ and insulatin^^H hosts, be it in the 
particle-hole symmetric limit!^ or away from it.l^ The 
LMA has so far been applied and adapted to a variety 
of impurity and lattice problems with strong electronic 
correlations, where its ability to cope with essentially all 
interaction regimes and to correctly describe all energy 
scales has proven very valuable. ^^ ■^^ In the present case 
of an impurity hosted in a narrow band, the LMA has 
to be generalized by using a renormalized, sum-rule com- 
pliant version of the original ladder-sum propagator for 
the transverse spin fluctuations. This extension is nec- 
essary, especially in the regime of moderate interaction 
strengths, to correctly capture the subtle low-energy dy- 
namics (see Sees. Ill A and IV B). 



The article is outlined as follows: after a brief presenta- 
tion of the model, Sec. [TT| defines the narrow-band regime 
and presents two different mean-field solutions. Sec. |III| 
starts with a review of the LMA, followed by a discus- 
sion of the physically relevant transverse spin fiuctuations 
and of two important sum rules for the associated po- 
larization propagator; its last paragraph is dedicated to 
the self-energy approximation implemented in practice. 



Sec. IV presents the LMA impurity spectra for different 
regimes defined by the strength of the on-site Coulomb 
interaction U; the spectra are compared to correspond- 
ing Hartree-Fock or perturbation theory results and, 
where available, to Hofstetter and Kehrein's Lanczos- 
determined spectra.!^ A Conclusion section closes the 
paper. 



II. HAMILTONIAN AND MEAN FIELD 
THEORIES 

The Hamiltonian for the AIM is given in standard no- 
tation by 

T U 

XI (^«kC+ Ck^ + h.c.) 
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where the first term describes electrons (of spin o ='f , ^) 
in a metallic host band of dispersion et. The following 
two terms refer to the impurity, with Si the impurity level 
and V the on-site Coulomb interaction. The final term 
describes the one-electron hybridization between the im- 
purity and host. 

Throughout this article, as in Refs. [2S] and 1301 the 
particle-hole symmetric AIM, obtained by setting Si = 



— [7/2, will be studied. In this case, the empty and dou- 
bly occupied impurity states are degenerate, whence for 
all interaction strengths the Fermi level remains fixed 
at its noninteracting value and the impurity charge is 
ni — {fii-f + fiij^) = 1. Regardless of the interaction 
strength, single-particle spectra are thus symmetric with 
respect to the Fermi level, w = 0. Persistent charge fluc- 
tuations guarantee the system's metallic character, allow- 
ing for the possibility to recast the exact single-particle 
impurity Green function as an infinite-order perturbation 
series, with each diagram depending solely on U and the 
noninteracting Green function 



g{uj) = [uj + iO+sgn(w) - A(tj)] ^ 



(2) 



In the latter expression, A(w) ~ AR(aj) — isgn(aj)Ai(w) 
stands for the hybridization function 



LJ + iO+sgn(a;) — 



(3) 



which condenses all relevant information about the host 
dispersion and the hybridization matrix elements Vn^. 

In this section, the single-particle impurity spectra 
of the AIM will be calculated in two different mean- 
field descriptions: Restricted Hartree-Fock (RHF), on 
the one hand, implements spin symmetry from the outset 
via identical impurity occupation numbers for both spin 
species; unrestricted Hartree Fock (UHF), on the other 
hand, seeks to determine the occupation numbers self- 
consistently — thus allowing for solutions with different 
impurity occupancies for f and 4,-spins at an intermedi- 
ate stage — and restores the full spin symmetry only at 
the very end. 

In both versions of the theory, the (causal) single- 
particle impurity Green function Giai'^) = ReC/io-(w) — 
i sgn(w)7rD°^(a;) can be deduced via standard techniques, 
e.g., the equation-of-motion method, after Hartree-Fock 
factorizing the 2-body term in the Hamiltonian ([I]), 
n,^% ~ ^^^t(%) + " yielding 



Giaiio) = + iO+sgn(a;) 



(4) 



where e^o- — £i + ^{n-i — cr/i) denotes the Hartree-Fock 
corrected impurity level. In the particle-hole symmetric 
AIM, where = —U/2 and rii = {fii^ + n.4 ) = 1, the 
Hartree-Fock corrected impurity level solely depends on 
the impurity moment fi = ( fii-^ — hi^ ), viz. 



U 



(5) 



with (T = + (— ) for t (i) spins. As a consequence of 
the Hamiltonian's invariance under spin inversion, solu- 
tions with a nonvanishing magnetic moment are doubly 
degenerate, — +\fJ,\ and — 

The present paper focuses on AIMs with narrow metal- 
lic host bands, whose width is much smaller than the hy- 
bridization strength at the Fermi level, Aq = Ai(a; = 0). 
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The physics of such a narrow-band AIM is naturally 
rather insensitive to the precise form of the hybridization, 
meaning that, without loss of generality, Ai(w) may be 
assumed to consist of a single flat band of intensity Ap, 
ranging from —D to with D <C Aq. As for any time- 
ordered Green function the real part follows via Hilbert 
transform, 



+ C30 

Ar(cj) = V I — 



duj' Ai(w') 



(6) 



so that in total 



In 



D 



D 



iAosgii{Lj)e{D-\uj\) . (7) 



A. Restricted Hartree-Fock (RHF) 

The spin-reversal invariance of the Anderson Hamilto- 
nian ([T]) implies that the average number of f and ^-spin 
electrons on the impurity has to be the same for any in- 
teraction strength. RHF theory acknowledges this fact 
from the outset by enforcing ( n,^ ) = ( n^^ ) for all inter- 
actions strengths, thus entailing = in eq. 0. As a 
result, RHF recovers the noninteracting Green function 
for both spin species and all interaction strengths: 



G^H^H EE 5H = [c. + iO+sgnH-AH]- 



(8) 



The associated single-particle spectrum D^^{uj) con- 
sists of two contributions: (i) a continuum for uj G 
[— £),£)], arising from the nonzero imaginary part of the 
hybridization; and (ii) two poles, one lying above and the 
other symmetrically below the continuum. The fraction 
of spectral weight residing in the poles depends strongly 
on the bandwidth of the host metal. In the usual wide- 
band model, D ^ Ao, most of the spectral intensity is 
concentrated in the single-particle band, while the poles 
are exponentially weak and hence irrelevant in practice. 

The situation is, however, radically different in the 
present narrow-band model, defined by _D <C Aq. Here, 
almost the entire spectral intensity resides in the poles, 
occurring at frequencies ±wo far outside the band. By 
analogy with an H2 molecule, these poles can be viewed 
as a bonding and an anti-bonding orbital, thus suggest- 
ing that the narrow host band to which the impu rity i s 
coupled behaves effectively as a single site or level.'^SESl 

In this case, the pole frequencies and weights can be 
obtained to good accuracy from eq. ([s]) by using the ex- 
pansion Ar(cli) ^ {2/tt)AoD/uj, valid for \uj\ ^ D: 



Wo 



^A„D 
8D 



(9a) 
(9b) 



The "bonding energy" ujq corresponds to the integrated 
hybridization, or the total hopping between host and 



impurityp^ subject to D ^ ^ Aq, it defines a second 
low-energy scale relevant in the narrow-band regime. 

In the following, it will sometimes be helpful to con- 
sider D and luq as independent parameters. In particular, 
in the limit of an infinitely narrow host band, 13 — >■ 0, this 
allows us to treat the host as a single level which couples 
via a finite bonding energy ujq to the impurity - a picture 
henceforth referred to as the two-site approximation. 

In addition to the orbital levels of the two-site approx- 
imation, the full noninteracting single-particle spectrum 
encompasses the aforementioned Fermi liquid continuum 
stemming from the hybridization band on the lowest en- 
ergy scale, |w| < D. Its integrated weight is of the order 
O (D/Ao) and thus weak in the narrow-band regime. 

The RHF description, naturally exact in the nonin- 
teracting limit, is expected to break down if U is much 
larger than the "molecular" bonding energy luq: in this 
case, the extra electron probed by is most likely to be 
introduced on an already singly occupied impurity which 
involves an energy cost of the order of the interaction 
strength U] the single-particle spectra will then be domi- 
nated by Hubbard poles separated by the Coulomb inter- 
action U rather than the molecular orbitals at a; = ±^9. 



B. Unrestricted Hartree-Fock (UHF) 

In UHF theory, the magnetic moment residing on the 
impurity will be determined self-consistently from the 
impurity Green function Qia- itself. For small interac- 
tion strengths, it is found to be zero, and UHF recovers 
the noninteracting solution g. This nonmagnetic solution 
becomes unstable above some critical interaction Uq — 
which turns out to be related to the "molecular" bond- 
ing energy ujq — and UHF then converges to a solution 
with a finite impurity moment. The Hamiltonian ([T]) is, 
however, still invariant under spin inversions, thus guar- 
anteeing for any mean-field ground state with positive 
moment /i = +/io the existence of another degenerate 
ground state with opposite moment, fi = —fiQ. Sub- 
sequently, quantum mechanical tunneling processes be- 
tween these mean-field ground states ensure their occur- 
rence with equal probability: 



G^"^(c^) 



1 



(10) 



Throughout the present paper and without loss of gener- 
ality, A- (B-) type impurities are assumed to be predomi- 
nantly t (4,) spin occupied, implying /io > 0. The spin 
symmetry of the Hamiltonian entails 5act(w) = Gb-(t{'-^), 
thus guaranteeing the spin independence of G^^^{uj) in 
eq. 

The necessity of the mixing process ( |10[ ) becomes par- 
ticularly obvious in the atomic limit, defined by vanishing 
hybridization matrix elements, T^ik = 0. In this limit, the 
impurity propagator is G^^{u}) — ^{[z+^]^^ + [z—^^]~^) 
[with z — Lu + i0^sgn(a;)]. Each of the two contribu- 



tions in eq. ( 10 ) yields one term of this exact result [via 
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A(w) = and = 1]: an f-spin electron can only be 
retrieved from an A- type impurity (first term), and can 
only be added to a B-type impurity (second term). 

In practical terms, the Green function t/A(T(w) — 
ReC/Acr('^) — isgn(w)7rZ?^^(w) is obtained from eqs. Q 
and ([5| with an impurity moment /i = +/io, calculated 
self-consistently from 



only retaining the pole contributions of D\^{uj) in the 
self-consistency equation (111. This yields the following 



expression which corresponds to the lower solid line in 
Fig.[l] 



{U^/Uf for [/ > C/o^ := 4a;o 

for u < m 



(13) 



MO = y du[D%{uj)-Dl^{^)\ ■ (11) 

— C30 

[^/B(t('^) follows equivalently for /i = — /io]- In Fig. [ij the 




FIG. 1: Magnetic impurity moment vs. Coulomb repulsion U 
for an AIM with bandwidth D = 0.01 Aq, where ujq ~ O.O8A0. 
Filled circles: numerically determined LMA moment ^, as re- 
quired by eq. ( 35 1; open circles: numerically determined UHF 
moment /ip. Solid lines: corresponding approximate analyti- 
cal expressions, eqs. (36 1 and (13 1. The critical interactions 
are U" 



O.23A0 for the LMA, and C/q 0.32Ao in UHF. 



self-consistent UHF moment /io, eq. (11), is plotted as a 



function of the interaction strength U for a narrow-band 
AIM with host bandwidth D = O.OlAo. 

The spectral density D^^(w) consists of a low-energy 
continuum for \oj\ < D, of net weight 0{D/ Aq), aris- 
ing from the finite imaginary part of the hybridization, 
and two pole contributions, one above the single-particle 
band, at = a;> > +D, and the other one below it, 
at cli = —i^< < -~D. Assuming both poles to occur far 
outside the single-particle band, the pole frequencies and 
weights can be obtained rather accurately by expansion 
of the hybridization function ([7|, AR(a;) ^ {2 / 1:) A^D / ljj: 



Wo 



l±2/o 



2/0 



2\/2/gTl 



(12a) 



(12b) 



with yo = ?7Aio/4a;o. 

For the above expansion to hold, ^ D is required, 
thus limiting the validity of eq. ( 12 ) for the low-energy 
pole at cj = — w< to interactions UfiQ <^ -Aq (while 
no restrictions follow for the high-energy pole w>). In 
this range of interactions, the renormalization effects 
in the single-particle continuum are still small and the 
UHF moment fj,Q can be obtained very accurately by 



For U < U§, the self-consistently determined moment 
vanishes and the UHF Green function coincides with the 
noninteracting or RHF solution, eq. ([S]). Above Uq, by 
contrast, a finite local moment forms on the impurity, 
and saturates rapidly as U is increased. With increas- 
ing /io, the pole at a; = cli> shifts rapidly away from 
its noninteracting value loq towards higher frequencies 
and gains in intensity. For U ^ U§, it becomes the 
upper Hubbard satellite at a;> ~ U/2, which overwhelm- 
ingly dominates the spectrum 15^^(0;) with a pole weight 

of (7> ~ 1 — O {[ujo/U]'^) . Simultaneously, the pole at 
Lu = —Lu^ moves from —loq towards the lower band edge 
—D and loses weight. For U§ U < Aq, its position 
and weight are given to good accuracy by 



2ul 



J 

U ~' 2 ' 

2ujoy 

u / ■ 



(14a) 
(14b) 



The first equation defines a third low-energy scale, J — 
Aujq/U, which lies between bandwidth and bonding en- 
ergy, D ^ J ^ wq, and accounts for the antiferromag- 
netic exchange between impurity and host. 

If U is increased to values of the order of Ao, the an- 
tiferromagnetic exchange J approaches the bandwidth 
D, and the above analysis — while still valid for the 
Hubbard level at w = uj^ — breaks down for the low- 
energy pole at uj — —Lo^. If, like in the present case, 
the hybridization band Ai(a;) has a discontinuity at the 
lower band edge —D, the logarithmic divergence of the 
related real part, AH(aj), still guarantees the existence of 
a low-energy pole; in the present model with a flat hy- 
bridization band it occurs, for U ^ Aq (where jig ~ 1), 
exponentially close to the band edge and carries expo- 
nentially small weight, which renders it insignificant in 
practice: 



D 



1 



2 exp 



ttU 
'2A^ 



2tiD 



q< 



■ exp 



ttU 
'2Aq 



(15a) 
(15b) 



Finally, according to eq. (10), the full UHF impu- 
rity spectrum can be obtained by superposing D\^{uj) 
and £'3^(0;), the latter of which follows by symmetry, 
D^^{uj) = D\^{uj). As illustrated in Fig. [Sj D^'^'^ {u) 
consist of a continuum for |cj| < Z), and two pairs of 
poles sX uj — ±a;> and lo — ±a;<. Before the onset 
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FIG. 2: UHF impurity spectra, TvAoD^^^^iuj) vs. w/Ao, for 
bandwidth D = O.OlAo, i.e., c^o ~ O.OSAo and Uq ~ 0.32Ao, 
and interaction strengths U/Ao = 0, 0.5, 1, 2 (from top to 
bottom). In the bottom panel, the Hubbard satellites at ~ 
±11/2 — ±Ao lie off plot range. Discrete levels are represented 
by vertical lines. The Fermi energy is cj = 0. 



curs qualitatively with second order perturbation the- 
ory in U (2PT),'32j and can be rationalized in terms of 
a simple two-site model in which the host is carica- 
tured by a single level coupled to the impurity.l^ De- 
spite these encouraging results, UHF suffers from several 
severe limitations: (i) in the two-site limit, 13 —> al- 
beit with finite ujq, exactly captured by 2PT, the low- 
energy single-particle levels are expected at w ~ ±3J/2 
for U ^ Uq, i.e., three times the corresponding UHF 
(ii) the UHF single-particle band, 
< Z3, is generally not a Fermi liquid: for 



pole frequency ( 14a 
present for \io 



interactions above U§, a nonzero impurity moment /ig 
forms, and the UHF zero-frequency behaviour D^^^{uj — 
0) = (l/7rAo)/[l + (C///o/2Ao)2] violates the Friedel sum 
rule^'^^' that compels the impurity spectra of the particle- 
hole symmetric AIM, for arbitrary interaction strengths, 
to be pinned at the Fermi level to their noninteracting 
value, D{uj = 0) = I/ttAq; and (iii), relatedly, due to the 
absence of dynamics UHF completely fails to capture any 
of the Kondo physics expected to govern the low-energy 
spectrum in the strong coupling regime, U 3> 4Ao. 

The present paper aims to transcend these shortcom- 
ings within the framework of the local moment approach 
(LMA). 



of moment formation, i.e., for U < Uq, the two pairs 
merge into the single pair of molecular orbitals shown in 
the first graph, at cj = iwg and of net weight q ^ 1/2 
each, and UHF coincides with the noninteracting solu- 
tion. As U exceeds Uq, the orbitals split up progressively 
into a stronger growing high- and a weaker growing low- 
energy component, as depicted in the second and third 
graphs of Fig. [2j In the second panel, where U = O.SAq 
is only moderately greater than U§ ~ 0.32Ao, the im- 
purity moment is already well established, /iq — 0.78, 
and the "Hubbard satellites" at cj ~ ±0.22Ao are not 
far from their terminal position, oj ~ zLU/2 — ±0.25Ao, 
carrying together more than 88% of the spectral inten- 
sity. Approximately a further 11% of the weight reside in 
the low-frequency poles at w ~ ±0.0 30Ao (expected at 
w ~ ± J/2 ~ ±0.025Ao from eq. (|14a[)), while the central 
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low-energy continuum carries only less than 1% of the 
spectral intensity. In the third graph, where U = Ag, 
roughly 98% of the spectral weight reside in the Hub- 
bard satellites at a; ~ ±0.49Ao, while the low-energy 
poles at w ~ ±0.015Ao appear close to the central con- 
tinuum, the latter being strongly renormalized from its 
noninteracting shape. Finally, these renormalization ef- 
fects are still enhanced for U = 2Aq (bottom panel), 
not yet fully in the strong coupling regime, resulting in 
a considerable violation of the Friedel sum rule (see be- 
low); here, almost all intensity resides in the Hubbard 
satellites (off plot range), and exponentially weak low- 
energy poles, with total weight 0.3%, are located slightly 
outside the central continuum [see eqs. (15)]. 

The abov e scenario , with a spectral evolution governed 
by ~ -\/2Ao-D/7r in weak coupling, and by U and 
J = Aujq/U in moderate coupling, C/q <C C/ < Aq, con- 



The LMAliSlMlssi expresses the impurity Green func- 
tion in a formalism employing two self-energies. At pure 
mean- field level this description reduces to UHF, dis- 
cussed in the previous section, with each of the self- 
energies arising from impurities with predominant t- or 
4,-spin occupation. For an impurity hosted by an insula- 
tor with a sufficiently large gap, this d oubly degenerate 
ground-state is actually observed.'2333Eil Conversely, for 
a metallic host, dynamical tunneling between the two 
broken-symmetry states can (and will) still lower the en- 
ergy and ultimately lead to the formation of a "Kondo" 
singlet ground-state wit h fully restored spin symmetry 
on the longest timescales! ^^ * ^^ ! The tunneling mechanism 
requires a (virtually) doubly occupied or empty impurity 
at some intermediate stage, implying that its rate ~ and 
the corresponding energy scale - diminish with increas- 
ing interaction strength. The LMA incorporates such a 
mechanism by coupling the single-particle dynamics to 
energetically low-lying flips of the impurity moment. In 
order to obtain a successful description of insulating and 
metallic phases within the same framework, this has to 
be done in a manner encompassing the possibility of self- 
consistent restoration of the spin symmetry at low en- 
ergies, as necessary for the preservation of Fermi-liquid 
behaviour on this scale. 

The implementation of the LMA follows Refs. |3T] and 
l32l For reasons analogous to those discussed in the con- 
text of eq. (10 1, the full impurity Green function G is 



again obtained by superposing A- and B-type impurity 
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propagators in a spin-rotationally invariant fashion 
1 



G(c.) 



[Ga.(w) + Gb.(w)] 



(16a) 



Making use of spin symmetry, GBa{^) = Ga-ct(w), yields 
the equivalent expression 

G{oj) = i [G„t(w) + G„^(^)] . (16b) 

Here, G is independent of the impurity type and the in- 
dices a may be suppressed for convenience, meaning that 
the individual (broken-symmetry) G„{uj) in eq. ( 16b I will 



henceforth be implicitly considered "A-type", i.e., of im- 
purity moment /i > 0, unless stated otherwise. Each 
of the Ggiuj) can be expressed in terms of a Dyson 
equation,'^ 



(17) 



thus defining the two self-energies Y,^{uj) and S^(a;) cen- 
tral to the present approach. Without loss of general- 
ity, the two self-energies can be separated into static and 
dynamic contributions, T,„{uj) — S^* + Y.^ri'-^)- diagram- 
matically, the former are suitably approximated by the 
Hartree tadpole (while the contribution of the Fock open 
oyster diagram vanishes), amounting to 



St* = 



u 



dc [D%[u:) - Dl^iu:)] 



(18) 



and the dynamical So.(w) are defined to contain "every- 
thing else" . Eq. (IT7|) may thus be rephrased as 



G^u) ^ [g-\Lo) - (St*-£.,) - I],(^)] ' . (19) 

This choice is particularly convenient if - as for the LMA 
- the dynamic self-energy contributions Ti„{uj) are to be 
diagrammatically constructed from the UHF propagators 
Q„ rather than from the noninteracting (7, since it entails 
the sum of the static contributions (in braces) to vanish 
i/the UHF moment /io, determined self-consistently from 
eq. (Ill, is used. 



The present two-self-energy description is by now well 
established and by no means exclusive to the LMA, but 
likewise emerges in other approaches as, e.g., in the NRG 
where it occurs for odd iterations of the renormaliza- 
tion group.f^Sl Nevertheless, it is desirable to connect it to 
the conventional single self-energy, defined by the Dyson 
equation G{<^) = [g-^{uj) - I](cj)]-i;3i 22 



S(^) 



Sf (cj) + S^(w) 





St(w) - S;(w) 


) 


2 







(20) 



Before specifying the class of diagrams to be retained 
to approximate the dynamical self-energies So-(w), the 
conditions necessary for Fermi-liquid behaviour to pre- 
vail at low frequencies shall be reviewed.'^ According 
to Luttinger,""^ the imaginary part of the single self- 
energy S(a;) is required to vanish as O (w^) at the Fermi 
level w = 0, which for the two self-energies So-(w) = 
S^(a;) — isgn(a;)Ej,(a;) employed by the present approach 
translates to the following two conditions: 



(^-0) 



O [<jj^) for w 



0; 



0) 



(21a) 
(21b) 



The fulfilment of the first condition shall be assumed 
from now on, and will be explicitly shown in Sec. |III 6] for 
the class of diagrams chosen in the following. The second 
condition requires the broken symmetry So-(w) to coin- 
cide with the fully symmetric single S(w) at the Fermi 
level, thus expressing the concept of a restored spin sym- 
metry on the lowest energy scaled which was alluded 
to in the beginning of this section. Here, the additional 
particlc-hole symmetry, S^(cj) = — S^(— w), in combina- 
tion with eqs. (21 1, require both self-energies S(t(w) to 
vanish at the Fermi level. This, in turn, automatically 
guarantees the Friedel sum-rule pinning of the spectra, 
D{ijj = 0) = I/ttAq. 

In its self-energy diagrams, the LMA incorporates 
transverse spin fluctuations responsible for a dynamical 
reversion of the impurity moment. At the simplest level, 
to be implemented in the following, this is achieved by 
coupling the ladder-sum polarization propagator to the 
single-particle Green function. 



where wavy lines represent the interaction U and left- 
going (right-going) solid lines stand for a UHF-like parti- 
cle (hole) propagator Q„. Note, however, that symmetry 
restoration is generally not automatically guaranteed by 
the above class of diagrams, but can be obtained through 
self-consistent determination of a free parameter. In the 
following, the impurity moment /i will serve as such, im- 
plying that it will generally differ at self-consistency from 
its UHF value /^o (see Sec. IIIB below). 

The above approximation for the self-energies 5]o-(w) 
has been motivated and discussed in detail in Ref. 1311 
I32I34I here, its basic properties shall only be briefly re- 
viewed: (i) in weak coupling, before the onset of moment 
formation, the UHF propagators coincide, for both cr, 
with the noninteracting Green functions g; up to (and in- 
cluding) second order in U, eq. (22) is hence equivalent to 
an ordinary diagrammatic perturbation expansion about 
the noninteracting limit, (ii) For arbitrary interaction 
strengths, the self-energy diagram (22) describes pro- 
cesses where after having added at i = 0, say, a a =4,-spin 
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electron to a — cr =t-spin occupied impurity, the original 
f-spin electron hops away, thus effectively flipping the 
impurity spin until its return at some later time t: this 
is precisely the dynamical spin-flip scattering, comprised 
in the second order of Anderson's poor man's scaling (see 
Ref. [9]), which allows us to restore spin symmetry and 
to capture the strong coupling Kondo limit, (iii) In the 
atomic limit, where l^k = and therefore Aq = 0, the 
self-energy (22 1 vanishes, since the — a-spin electron can- 



the right-hand side (rhs) of eq. ( 24 1 , yield^^ 



not leave the impurity at all. The LMA reduces to simple 
UHF, and thus becomes exact in this limit, predicting an 
impurity moment of /iq = 1 and single-particle poles at 
uj = ±U/2. This is salutary, since the atomic limit may 
be considered as an extre me of the local moment phases 
observed in a gapped AIM}22M2Hlllit also plays an impor- 
tant role in the context of the breakdown of the skeleton 



expansiorF^ taking place already in the metall ic regime 
close to a Mott metal-insulator transition.'^SMl 

Before discussing the LMA impurity spectra in Sec. |IV[ 
the spin fluctuations and their dynamical coupling to the 
single-particle Green function will be examined in the 
following two paragraphs. 



A. Transverse spin fluctuations 

As already mentioned in the previous paragraph, one 
of the main ingredients for the LMA self-energy is given 
by dynamical spin-flip scattering processes opening the 
possibility to restore the spin symmetry on the lowest 
energy scale. Embodied in the transverse impurity-spin 
polarization propagators, Tl+-{t) = i ^ f {S'+(t)S'r(0)} ^ 

and n ^{t), these scattering processes are at the sim- 
plest level accounted for by a sum of repeated particle- 
hole interactions, as depicted by the ladder bubble in the 
rightmost self-energy diagram in eq. (22 ) which, for a 



{—a =t), leads to the familiar expression 



n+-(c.) = 



"n+-(^) 
1 -c/on+-(w) ' 



(23) 



where 



+00 



°II+-{lu) = i 



duj' 



g^{j)g^(uj' ~uj) (24) 



-Im°n+-(a;) = 6'(w) \ duj' DI^uj') D^{uj' - uj) 



+0{-Uj) J duj' Dl{uj')D'l{u;' -Lo). 
-H 

(25) 

The first term describes a spin-flipping particle-hole ex- 
citation where an electron is removed from an occupied 
t-spin state below the Fermi surface (w' — w < 0) and 
placed into an unoccupied 4,-spin state above the Fermi 
surface (w' > 0). Similarly, the second term accounts for 
the recombination of an f-spin electron into an empty 
^-spin level below the Fermi surface. 

In the present AIM with a narrow host band, these pro- 
cesses operate mainly between the UHF orbitals, giving 
rise to poles in ;^Im°n'' {uj) at w = 2a;> and uj = — 2a;< 
with respective pole weight . 

For a vanishing impurity moment, as appropriate for 
weak coupling, the UHF levels coincide with the "molec- 
ular" orbitals of the noninteracting limit, entailing poles 
in ilm'^n+~(a;) at frequencies w = ±2ojo and of spec- 
tral weight q ~ 1/4 each. Conversely, in moderate to 
strong coupling, particle-hole excitations between the 
lower and the upper UHF Hubbard level become domi- 
nant, producing a high-frequency pole at = 2a;> ~ [/ of 
weight g?. ~ 1 in -Im'^H^ (cj), and a corresponding low- 
frequency pole at Lj — — 2aj< whose intensity decreases 
rapidly (first algebraically for If^ <^ U < Aq, where 
^ [2ujo/U]^, and then exponentially for U ^ 4Ao). 

In addition to the pole contributions -Im°H^ (w) 
contains three weak continua stemming from processes 
involving the narrow UHF single-particle band at the 
Fermi level. Two of these continua have their origin in 
electronic transfer between the single-particle band and 
one of the UHF levels; they are found at uj £ [w> , a;> +D] 
and at a; S [— w< — D, — w<], with respective net weight 
O (g^D/Ao), and turn out to be of httle importance in 
the following. The third continuum appears at low fre- 
quencies, \u!\ < 2D, and owes its existence to particle-hole 
excitations within the narrow UHF single-particle band 
itself. According to eq. (25 1, its behaviour in the vicinity 



of the Fermi level is essentially governed by the overlap 
of the two D° (a;)-bands; for a metallic host, the intensity 
of the latter being finite at the Fermi level uj = 0, this 
results in a continuum vanishing linearly for small w, 



is the bare transverse-spin bubble. [The second po- 
larization propagator follows by particle-hole symmetry, 

n+-(w) = H-+(-w).] 

Using the Hilbert transform relation (|6| for Gaiuj) on 



ilm°H+-(c.) \uj\Dl{0)D',{0) , (26) 



a prop erty that is needed to prove the fulfilment of con- 
dition (21a) by the present theory.ISS 
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1. Stability of the ladder sum propagator. 



2. Sum rules and renormalized ladder-sum propagator. 



From the Hilbert transform relation ([6| for (cj), 
and the positivity of its imaginary part, one deduces an 
analyticity condition for the polarization propagatorpS 



+ 00 

Ren+-{uj = 0) = / — 



duj' Imn+-(cj' 



> , (27) 



which naturally holds for and the exact 11+ . 

Conversely, for an approximate polari zati on propaga- 
tor, like the present ladder-sum , eq. ( p?] ) is not auto- 
matically satisfied, but rather constrains the applicability 
of the approximation to a certain range of parameters U 
and /i. 



According to eq. (23) and granted Im°n^ = 0) = 
0, the behaviour of Rell^ (w = 0) is solely controlled 
by the bare Re^II'' (a; = 0). For a vanishing impurity 
moment, as enforced in RHF and found self-consistently 
in UHF for U < = 4wo, Re°n+-(a; = 0) ~ l/4wo 
follows to good accuracy from the two-site approxima- 
tion, implying that the corresponding ladder-sum \V 
satisfies the positivity condition ( 27 ) . 



Above C/q, the /i = ladder-sum TV — which cor- 
responds to an ordinary random phase approximation 
(RPA) in the transverse spin channel — becomes un- 
stable to spin excitations since the underlying mean-field 
theory predicts a transition to a phase with a finite mag- 
netic moment. A finite /i, on the other hand, all ows for 
a direct calculation of Re^II^ {uj = 0) from eq. (24) by 
means of the identity Q^{ijj) — G^ii^) = —U^Q^{uj)Q\{uj) 
and of particle-hole symmetry = D° ^(— w): 



i7Re°n+- 



0) 



At 



The latter equation shows that if the moment is de- 
termined self-consistently from the UHF propagators, 
fi = Ho as given in UHF by eq. (11), the rhs evaluates to 
unity and the corresponding ladder-sum H"' , eq. ( 23 ) , 
has a pole at a; = 0. While appropriate for an insulator, 
where flipping the impurity spin costs no energy, in a 
metallic host such a spin fiip is expected to be governed 
by the Kondo effect and thus to involve a finite energy. 

In the framework of the LMA, this picture emerges nat- 
urally if /i is increased above its UHF value fiQ, thereby 
shifting the pole in H^ (w) from a; = to a small but 
positive frequency lo^ — closely related to the Kondo 
energy — and rendering H'' analytic in the sense of 
eq. (p7). 



Any acceptable approximation for the spin-flip polar- 
ization propagator should fulfil the following sum rules: 



+00 



± 



dcj 



ImH+ (w) = ( c+Ci^,c+Ci^ 



" V 

(29) 

which naturally hold for the noninteracting ''H"' (cj), 
and likewise for the exact H"' (w). 

The second sum rule reflects the fact that, in weak 
coupling, ImH"' (w) is symmetric about the Fermi level 
since the impurity is found to be occupied by electrons of 
both spin species with equal probability; after the onset 
of moment formation (and before symmetry restoration) , 
the impurity is predominantly f-spin occupied and spec- 
tral intensity at positive energies, associated with pro- 
cesses flipping the impurity spin from f to J,, should dom- 
inate. 

As for the first sum rule, its rhs is comprised between 
1/2 and 1 for the particle-hole symmetric AIM, since the 
expectation to find a doubly occupied impurity ( n^^n^ ) 
varies from 1/4 in weak to in strong coupling and, in 
any case, n^^ + rii^ — 1; both limits are captured by 
the approximate expression |[1 + jJ?] which follows by 
recasting the double occupancy as ( ni^hi^ ) = j{{fii} — 
{fj'i}) and subsequently factorizing the quartic Fermion 
operators via a Hartree-Fock decoupling. 

For the present narrow-band AIM, the above sum rules 
are naturally fulfilled in weak coupling, where H^ (w) 
reduces essentially to the noninteracting "H^ (w). More 
surprisingly, also in the strong coupling (Kondo) regime, 
U ^ 4Ao, both sum rules are found to be approximately 
satisfied for values of the magnetic moment required by 
symmetry restoration. 

Problems are found to arise primarily for moderate 
interactions, large enough to find the impurity moment 
well established, but far from the Kondo limit. For this 
range of interactions, the main features of the ladder-sum 
polarization propagator H^ (w) are accurately described 
in the two-site approximation; it reduces to two poles, 
with frequencies and weights given analytically by 



= 2wn 



iy±\/hy)^ + i 



7y ± Vi^y)^ + 7 



where y = fj,U /Uq and 7 



1] [hy? 
= 1 - 



+ 7] 



(30a) 



(30b) 



Inciden- 



tally, eqs. (30) encompass the bare (broken symmetry) 
°H'' (cj), which follows for 7 = 1, as well as the nonin- 
teracting polarization bubble, obtained by setting U — 
and 7 = 1. 
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Within the two-site approximation, appreciable sum- 
rule violations occur for U ^ Uq, where the stability 
criterion (27 1 and the definition of the magnetic moment 



itself constrain /i to the narrow interval between the UHF 
moment fj,Q and 1. In this regime, the deviation of the 
moment from its UHF value defines a small positive pa- 
rameter S — — /Zq. The position and weights of the 
ladder-sum poles, eqs. (30), can then be expanded in 



terms of the square root of this parameter which, using 
7 ~ /io5 + 0{6^), yields 



(31a) 
(31b) 



By consequence, the integrals in eq. (29) evaluate to + 
'0{V6) and q+~q^ = /xq; the l/y/S- 



9m = 2a;o/C/VA*o^ 
singularity of the first result indicates a strong violation 
of the first sum rule (while the second is found to be 
fulfilled to leading order even in this regime). 

Within the present approach, this problem can be over- 
come by renormalization of the ladder-sum propagator, 
i.e., by multiplication of Imn^* (w) with two different 
constants above and below the Fermi level, w = 0, cho- 
sen to comply with both sum rules (29). [The corre- 



sponding real part can subsequently be calculated from 
the Hilbert transform In the two-site approxima- 

tion, the renormalized weights — which have to be used 
instead of eq. ( 30b ) — are explicitly given by 



(32) 



Henceforth, unless explicitly stated otherwise, H"' (w) 
will always stand for such a renormalized spin-flip ladder 
sum. 



3. The full spm-flip polarization propagator. 

Beyond the two-site approximation, where the polar- 
ization propagator solely consists of two poles, with fre- 



quencies given by eq. (30a) and renormalized weights 
( 32 ) , the full ladder-sum H^ contains the following ad- 



ditional features: 

(i) The full H^ (oj) inherits its continua from the 
bare °U'^^{uj); these are located at frequencies cj e 
±[u}^,uj^ +D] and |w| < 2D. 

(ii) For moments fi which only exceed the correspond- 
ing UHF moment very slightly, as appropriate in strong 
coupling, U >• 4Ao, the pole frequencies (30a) may he 
within the low-energy continuum |a;| < 2D; instead of 
poles, the full H^ (ui) will thus have sharp resonances at 
these frequencies. 

(iii) In addition to the above mentioned poles, the full 
n+^(w) shows a third collective pole, induced by the 



RPA-like structure of eq. ( 23 ) and the logarithmic singu- 
larity of Re°H^ (uj) at the upper edge of the polariza- 
tion band found for uj g [— Ci;< — D, — w<]; of tiny weight 



throughout the whole range of interactions, this pole can 
be safely neglected in practice. 



B. The LMA self-energy 



For a 4,-spin electron, the self-energy diagram (22) 
translates to 



+ 00 



d( / 

S;(^) = / — H+-(c.')t/T(^-^') ; (33) 



and the f-spin self-energy follows by particle-hole sym- 
metry S^(w) = -E^(-a;). 

In analogy to the calculus of the spin-flip polariza- 
tion propagator, useful expressions for the self-energy 
I]^(a;) = I]^(w) — isgn(a;)E|(a;) are obtained by insert- 
ing the Hilbert transforms for the broken-symmetry H'^ 
and in the integrands on the rhs of eq. (33), yielding 
for the imaginary part 



S]_(w) = 0(w)t/2 / duj'Imn+- {uj')D^{uj-uj') 



+0{-uj)U^ J duj'linn+-{uj')D^uj ~Lu'); (34) 



and a corresponding real part, I]j^(ti;) 
Hilbert transform. 



which follows by 



The evaluation of the integrals in eq. ( 34 ) is greatly fa- 



cilitated by the structure of the polarization propagator 
H"' (w), whose spectral intensity resides — as pointed 
out in Sec. 



throughout the entire range of inter- 
actions almost exclusively in two sharp low-energy modes 
at a; = uj^. Except in strong couphng, U 4Ao, the 
dominant contribution to the self-energy arises by cou- 
pling these modes to the UHF f-spin orbitals, resulting 
in two poles or sharp resonances in S|(a;), the first oc- 
curring at w = w+ + iuj< > (net weight U'^q^q^), and 
the second at a; = uj^ — w> < (net weight U'^q~qy). 

The self-energy pole at positive frequencies is associ- 
ated with single-particle excitations where an extra J,- 
spin electron is introduced on the impurity, and the initial 
f-spin impurity electron hops off to the host. The same 
result could have been obtained from the original ground- 
state by placing the extra electron directly in the empty 
t-spin UHF orbital (energy cost a;<) and simultaneously 
flipping the impurity moment from f to 4^ (energy cost 
u!^). Analogously, the E|-pole at negative frequencies is 
caused by a J,-spin electron being taken off the impurity 
and subsequently replaced by an f-spin electron from the 
host, yielding an f-spin electron in the UHF orbital be- 
low the Fermi level (energy gain w>) and a flip of the 
impurity moment from f to f (energy gain As this 

requires the impurity to be initially f-spin occupied, the 
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latter self-energy contribution will decline as the impu- 
rity moment fj, approaches saturation. 



In addition to these collective modes, eq. (34) predicts 



several sets of narrow self-energy continua that arise from 
band contributions in Z55 and 11^ : such continua occur 



for Lo > o;,^, for lu < 



and for \uj\ > a;> -I- a;< (all of 
width D), for \uj\ > uj^ (width 2D), and for \uj\ < 3D. 
All these continua remain generally weak over the whole 
range of parameters, albeit with one exception: as the im- 
purity moment /i approaches its UHF value from above, 
H Ho + 0, and consequently ui^ — )■ — a situation 
that becomes relevant in the strong coupling Kondo limit 
discussed in Sec. |III B 2| below — the minor UHF single- 
particle pole vanishes exponentially [see eq. (15b)], and 
the dominant self-energy contribution is transferred to 
the low-energy continuum located at |a;| < SD. 

That the low-frequency behaviour of the latter con- 
tinuum fulfils the first condition for symmetry restora- 



tion, eq. (21a), may be seen as follows: the continuum is 



generated, as in eq. (34), by convoluting the UHF band 



for f-spin electrons with the low-frequency continuum of 
IniH"' ; the latter vanishes linearly in a; as w — >■ [a 
consequence of eq. ( 26 ) in combination with the analyt- 



icity condition (27), which for a ladder-sum propagator 
implies 1 - f/Re^+"(a; = 0) > 0], entaihng the ir con- 
volution to behave as [S|(w) =]S]j^(ti;) cx for \uj\ <C D. 
The second condition, eq. (21b), by contrast, is not 



automatically fulfilled by the present theory. As pointed 
out before, it reduces under particle-hole symmetry to 
tf{u} = 0) = ^* + T.f{uj) = 0, which, under the as- 
sumption Ei'*^* ~ Ufi/2, becomes 



Ef{u = 0) + 1^ = 



Within the LMA, this condition will be satisfied by 
tuning the impurity moment fi within the interval fiQ < 
< Ij prescribed by the analyticity condition (27). De- 



spite the generally narrow range of possible values, a so- 
lution can always be found since the dynamical contri- 
bution, Ej^(w = 0), is very sensitive to the exact posi- 
tion of the majority self-energy pole and hence implic- 
itly to the spin- flip scale w^. The latter quantity is 
highly responsive to changes of the impurity moment 
and thus effectively controls the condition ( 35 ) . In a 
more general perspective, o;^ — which vanishes for an 
ordinary UHF ground-state with finite impurity moment 
(see Sec. HI A 1 ) — may be considered as a new order pa- 
rameter that, for a state with local moments, determines 
whether Kondo physics takes place or not. 

The above analysis is corroborated by the numerical 
results displayed in Fig. [l] confirming that, for a wide 
range of interactions, U ^ Uq, the impurity moment fi 
needed to comply with condition ( 35 ) is indeed very close 
to the UHF moment ^o- 

In the remaining two paragraphs of this section, condi- 
tion ( 35 ) in combination with the LMA self-energy will be 



worked out: (i) the two-site limit, appropriate for weak 
to moderate interactions; and (ii) the strong coupling 
Kondo limit. 



1. The self-energy m the two-site limit. 

The two-site approximation — Z? — >■ albeit with finite 
Wo ~ ■y/2Ao£'/7r — allows for the LMA self-energy to be 
written in simple analytical terms and, more importantly, 
suggests a two-pole structure for Y,a(uj) that, for w S> 
\D\ and up to moderate interactions U < Aq, correctly 
reproduces the main results of a complete version of the 
LMA for the present narrow-band AIM. 

In weak coupling, [/ <C J/g, UHF predicts a solution 
with a vanishing impurity moment, /xg — 0, and single- 
particle propagators which, independent of spin and im- 
purity type, reduce to the noninteracting Green function, 
Qa{uj) = g{uj). The corresponding self-energies T,„{uj) 
are hence odd functions of w and their real parts sat- 



isfy the condition ( 35 ) by symmetry. As obvious from 



eq. (22 1, the leading contribution to the self-energy dia- 



grams stems from ordinary 2PT about the noninteracting 
ground state and is hence overwhelmingly dominated by 
poles at w ~ ±3wo of weight Q ^ U'^/8 each. 

Above some critical interaction If^ — which is of the 
same order of magnitude but slightly smaller than the 
corresponding UHF critical interaction U§ — the non- 
magnetic solution becomes unstable and a finite mo- 
ment forms on the impurity. As illustrated in Fig. [T] 
the LMA moment /z saturates rapidly with increasing U, 
and its numerical values fit accurately to a square-root 
law [which, in contrast to and despite the similarity with 



(35) eq. (I3|, is not even exact in the two-site limit]: 



y^l - {U'^/Uf for U >U'' := 2^2 



UJQ 







for U <U'' 



(36) 



As pointed out above, symmetry restoration and con- 
sequently the Fermi-liquid nature of the single-particle 
excitations hinge, in this regime, almost exclusively on 
the energy cost Wm for flipping the impurity moment. 
For simplicity, in what follows, the latter quantity will 
be thought of as an independent parameter, to be de- 
termined from condition ( 35 ) , and the impurity moment 



will be kept at its UHF value fio instead. 

Under these assumptions (which are in conc ord to 
leading-order with the analysis in Sec. Ill A 2), H^ 
is dominated by poles at w,^ = ii^m, of renormal- 
ized weights = ^{1 ± ^o)^) and E^(a;) is consti- 



tuted by a majority pole at oj 



+ w< (net weight 



U^q^q^ = 4wQq>) and a minority pole at a; = — 0;^ — 
(net weight U^q^q> = 4a;oq<). The magnetic energy 



scale follows finally from eq. (35), viz 



studied in two regimes for which analytic results can be 




(37) 
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where, again, J = Aluq/U stands for the antiferromag- 
netic exchange coupling constant. 

At the critical interaction for moment formation in 
UHF, where U = Uq and U/A = J = ujq, the magnetic 
energy starts out at = V^ujq, and then decreases as 
uJm ^ 3J/2 for interactions U ^ Uq (the relevance of 
the two-site approximation being subject to lo^ ^ D, or 
equivalently U < Aq as previously). 

Due to the simplicity of the two-site approxima- 
tion, two further virtues of the LMA become mani- 
fest: first, the renormalization procedure of IT"' (uj) 
— which enforces the sum rules ( 29 1 — guarantees 



the LMA self-energy ( 33 1 to fulfil an analogous set 



of sum rules for the interaction self-energy 



+00 
± 





da; ~T , , 



(38) 



the sole condition being that the many-body expectation 
values on the rhs are evaluated by a Hartree-Fock factor- 
ization, as shown on the second line of eq.(38) (which, 
again, correctly captures the limits of weak and strong 
coupling) . 

The second virtue concerns the conventional single self- 
energy E: in Lanczos calculations by Hofstetter and 
Kehrein,'2Sl this E(a;) is found to have poles on an energy 
scale UJ ~ y/ AqD which can not be explained in any or- 
der of the skeleton expansion.l^^ISfil Moreover, these poles 
have been shown to depend little or not at all on the 
interaction strength,'351 occurring in the limits of strong 
and weak coupling exactly at w = ±3^0, with net spec- 
tral weight Q ~ each. Responsible in the atomic 
limit t^ik ^ LdQ = for the E(a;) ^ l/ui characteristics of 
the insulator, they are an intrinsic feature of the narrow- 
band AIM0 

It has been argued in the beginning of this paragraph 
that, in weak coupling {U < If^), both LMA self-energies 
coincide and inherit their main properties from 2PT, 
whose poles in turn precisely respect the desired prop- 
erties. But also for U ^ , where the impurity moment 



saturates and w,, 



-,J, such poles occur on a similar 



energy scale in the symmetry-restored single LMA self- 
energy, E(w). The latter depends, as in eq. (20), on the 
self-energies So- which in the limit U ^ reduce to 



UJ \uj + 



AujI 



(39) 



and Ti-^iuj) = — E^(— w); insertion of which into eq. (20) 



generates a single self-energy constituted, to leading or- 
der in J '-^ 1/C^, by poles at a; ~ ±y/bujQ of spectral 



weight Q ~ ?7^/8 each. Hence apart from the pole fre- 
quency prefactor (-^5 ~ 2.23 instead of 3), the LMA 
self-energy matches the above mentioned properties of 
the exact solution also for U ^ W^. The discrepancy in 
the prefactor is possibly due to the inability of the two- 
site approximation to describe the true strong coupling 
regime, U » 4Ao. 



2. The self-energy in the Kondo limit. 

In genuinely strong coupling, U ^ 4Ao, the polar- 
ization propagator 11"' (uj) consists in essence of a sin- 
gle sharp resonance peaked at a frequency uj^ which is 
tiny compared to any other frequency scale involved in 
the problem, even D. The corresponding self-energy, 
obtained by convoluting this resonance with the single- 
particle UHF propagator as in eq. (34), is thus governed 



by resonant spin-fiip scattering within the metallic single- 
particle band at the Fermi level itself, i.e., the Kondo 
effect, whereas the orbital contributions, predominant in 
the two-site limit, vanish exponentially [see eq. (15b)]. 
Under these circumstances, and granted /ito — 1, 



condition (35) can again be solved analytically, yield- 
ing an exponentially small spin-flip energy, uj^ ~ 
D exp[— 7rC//8Ao], which — apart from the prefactor — 
concurs with the LMA results for the AIM with a flat 
and infinite wide hybridization band.l^Far from the band 
edges of the low-energy continuum, the interaction self- 
energy can be written as a function of a single parameter, 
UJ — UJ /uj-i-a, viz. 



i;^(w) 4Ao 



-\n\l-Oj\ -i6{Cj-l] 

TT 



(40) 



and is otherwise independent of the original parameters 
U and D. 

This low-fequency scaling behaviour of the self-energy 
- which is ultimately responsible for the scaling prop- 
erties of the single-particle resonan ce - is a hallmark 
of AIMs with metallic hosts. 'i^lMiH xhe associated pre- 
diction of an exponentially small magnetic energy scale, 
closely related to the Kondo temperature, emerges sim- 
ilarly from Anderson's poor man's scalin^^ for the s-d 
model, and from a strong-coupling expansion of the ex- 
act Bethe ansatz solution for the AIMM 



IV. SPECTRAL EVOLUTION IN THE 
NARROW-BAND REGIME. 

In this section, the numerically obtained LMA im- 
purity spectra D{uj) ~ — ^sgn{uj)lm G (uj) [with G{uj) 
the spin-symmetric single-particle Green function of 
eq. (16b)] will be presented next to corresponding UHF 
and 2PT results. The section also comprises a compar- 
ison of the LMA spectra with Lanczos calculations per- 
formed by Hofstetter and Kehrein^^ for an 11 + 1-site 
Anderson star with host bandwidth D = IO^^'Aq. 
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Further analytic rationales reveal how the spectrum 
evolves with interaction strength and suggest an interpre- 
tation of its main features in simple physical terms: for 
weak to moderate interactions, U < Aq, the reasoning is 
based on the two-site approximation where the host band 
is taken to be infinitely narrow; conversely, in strong cou- 
pling, U ^ 4Ao, similarities with the opposite case of an 
infinitely wide host band emerge in the immediate vicin- 
ity of the Fermi level, and ultimately lead to the charac- 
teristic scaling behaviour of the Kondo resonance!^^'^ 



A. Weak coupling: U ^ W 

For U much smaller than the critical interaction for 
moment formation, W^, the LMA connects smoothly to 
2PT (see Sec. IIIB 1 ), which in turn exactly captures the 
two-site limit for any interaction strength. By analogy 



to the latter, extensively discussed in Ref. [301 the LMA 
single-particle spectrum (depicted in the last graph of 
Fig. [3]) is overwhelmingly dominated by the "molecular 
orbitals," occurring at minimally lower frequencies |w| < 
Wo than in the noninteracting limit. Additionally, for 
any nonzero U, a pair of weak poles arises at cj ~ zbSwo: 
these are the precursors of the Hubbard satellites, and 
start out with net spectral weig ht q^O {[U/uJo]^) each. 
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FIG. 3: UHF (top), 2PT (middle), and LMA (bottom) impu- 
rity spectra, tvAqD^uj) vs. uj/Aq, for bandwidth D = O.OlAo, 
i.e., cjo — O.O8A0 and If^ ~ 0.23Ao, and interaction strength 
U = O.IAq. Discrete levels are represented by vertical lines. 
The Fermi energy is = 0. 



(middle panel), namely: (i) the "molecular orbitals", oc- 
curring at w ~ ±0.077Ao, dominate the spectrum and 
carry in total more than 92% of the intensity; (ii) the 
outer set of poles, at w ~ ±0.22Ao, carrying approxi- 
mately 6.6% of the total spectral weight, occur slightly 
inside their 2PT counterparts at w ~ ±0.25Ao; (iii) vari- 
ous weak band contributions which (with one exception) 
also bear great similarity with 2PT, the most prominent 
being the essentially unrenormalized Fermi-liquid con- 
tinuum of width 2D and net spectral weight O (D/Aq) 
around the Fermi level. These results clearly contrast 
with UHF, in the top panel of Fig. [3] whose spectra, for 
the reasons already pointed out in Sec. |HB coincide with 
the noninteracting limit. 



B. Moderate coupling: U" <. U < Ao 

In the regime of moderate interaction strengths — de- 
fined hy Lf^ ^ U < Ao or, equivalently, D <^ J <^ ujq — 
a well-established local moment resides on the impurity. 
The energy cost Wm — for hipping this moment being 
large compared to the host bandwidth D, the two-site 
scenario of Sec. |IIIB 1| remains appropriate — except, 
of course, within the low-energy Fcrmi-liquid continuum 
\uj\<D itself. 
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FIG. 4: UHF (top), 2PT (middle), and LMA (bottom) impu- 
rity spectra, tiAqD^u) vs. u/ Aq, for bandwidth D = O.OIAq, 
i.e., cjo — O.O8A0 and ~ 0.23Ao, and interaction strength 
U = O.5A0 (J — O.O5A0). Discrete levels are represented by 
vertical lines. The Fermi energy is a; = 0. 



Fig. [3] depicts the impurity spectra for D — O.OlAo 
and U — O.IAq, i.e., an interaction strength that is weak 
compared to Aq, but already appreciable with respect 
to U" ~ 2^/2uJo ~ 0.23Ao (or = Aluq ~ 0.32Ao). 
Since U < W^, the LMA converges to a solution without 
magnetic moment, and the numerically computed LMA 
spectra (bottom panel) carry mainly the signature of 2PT 



For a host of bandwidth D ~ O.OIAq, chosen in Fig. |4] 
for clear resolution on all relevant energy scales, the 
above conditions of moderate interaction strengths are 
never truly met, since ~ 0.23Aq and Aq are of sim- 
ilar orders of magnitude. Nevertheless, for U = 0.5Aq 
(as in the figure), where the LMA impurity moment is 
well established (fj. ~ 0.89 vs. /iQ ~ 0.78 for UHF, 
see Fig. II]), the moderate coupling regime of the two- 
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site approximation is a good starting point for the anal- 
ysis of the impurity spectra. Their most prominent 
feature, accounting for approximately 68% of the spec- 
tral intensity, are the high-energy Hubbard satellites at 
~ -I- I J ~ O.33A0. They are found to occur re- 
markably close to their 2PT counterparts at |w| ~ 0.35Ao 
(expected in the two-site limit of 2PT at |a;| ~ 
0.38Ao), whereas the UHF correction to the atomic limit 
(where \uj\ = operates in the opposite direction, 
~ — i J ~ O.22A0. A new feature on the high en- 
ergy scale are the side bands for |w| > a;> -|-a;< ~ 0.27Ao 
(and, very faint, for \uj\ > lu^ ~ 0.25Ao), stemming 
from many-body self-energy continua (see Sec. IIIB): 
strongly enhanced with respect to their 2PT progenitors 
(located in the middle panel at = 2wo — O.I6A0 and 
Wo — O.O8A0), these side bands will, for larger U, eventu- 
ally merge with the Hubbard levels, thereby destroying 
their discrete nature. This, in turn, may be viewed as 
an extreme version of the many-body broadening effect 
generic to the AIM a nd likewise observed for hosts of 
large band width. l^iEIl 

In the low-energy sector — governed by the antiferro- 
magnetic scale J = 4u!q/U ~ 0.05Ao and the bandwidth 
D = O.OlAo — the LMA impurity spectra reveal a much 
richer structure than their UHF or 2PT counterparts. 
While all three approaches predict the low-energy con- 
tinuum, ranging from uj = —D to D and accounting for 
0.4% of the total spectral intensity, to suffer very little 
renormalization with respect to the noninteracting limit, 
the spectra differ significantly on the J-scale, where sim- 
pler UHF and 2PT both predict a single pair of discrete 
levels (which can be viewed as remnants of the molecular 
orbitals) . An additional, second pair of levels emerges in 
the LMA already in the framework of the two-site ap- 
proximation: the pole equation for the 4,-spin electron, 
Lj — Ar (w) — S^(w) = 0, with S^(w) from eq. 
the upper Hubbard satellite together with two 
poles at, to leading order in J, 



39 ) , yields 



ow-energy 



1± ^/l7 



J 



(41) 



while, the ^-spin pole equation, leads to the lower Hub- 
bard satellite and poles at w = — cji.2. 

The low-energy poles at w = ±wi^2 can be interpreted 
in the following way: the outer levels at |cj| = \lui\ ~ 
1.28J ~ 0.065Ao (of net weight qi ~ 0.51[J/wo]2), re- 
sult from a shift of the UHF orbital remnants, initially 
found at \uj\ = aj< — 5 J, and occur remarkably close to 
their 2PT counterparts at |cl;| ~ | J ~ O.OTGAq; the inner 
poles, at frequencies |w| = — 0.78J ~ 0.04Ao and of 
net weight (72 — 0.12[J/wo]^, are entirely new features, 
which owe their existence to resonant spin-flip scatter- 
ing as accounted for by the majority self-energy pole at 
Lu = a;< + ujni — 2 J. In UHF, this pole is naturally ab- 
sent due to the total lack of dynamics; but it also misses 
in 2PT, albeit for subtler reasons: assuming a perfect 
parity between host and impurity, as appropriate for an 
actual two-site problem, a resonance process that tends 



to localize a spin flip on the impurity rather than on the 
host is ruled out in 2PT by construction. 

Despite the rather poor concord with the definition 
of moderate coupling, the numerically calculated LMA 
spectrum of Fig. |4] corroborates the above analysis to 
reasonable accuracy, displaying in the low-energy sector 
a four-peak structure with (i) levels at ~ O.OSAq ~ J 
(instead of the expected 1.28J) and carrying about 21% 
of the spectral intensity; and (ii) sharply peaked continua 
at \uj\ ~ 0.7J ~ O.O35A0, accounting for roughly 8% of 
the spectral weight, stemming from the merger of the 
above mentioned inner poles, expected at a; ~ ±0.78J, 
with self-energy bands. 



C. Comparison with Lanczos spectra. 

In their article,'251 Hofstetter and Kehrein present low- 
frequency Lanczos spectra for an Anderson impurity cou- 
pled to a host of 11 sites, with a bandwidth of 13 = 
10^^ Ao, entaihnga "bonding energy" ofwo ~ O.OO8A0 = 
8QD. They study two different interaction strengths: the 
first, U ~ O.2A0, implies J ~ 12.7D much larger than 
D but much smaller than If^ ~ 230D, thus matching 
the above definition of moderate coupling; the second, 
U — 4Ao, implying J ~ 0.6313 slightly smaller than the 
bandwidth, is in the crossover region between moderate 
coupling and the Kondo regime. Fig. [5] shows LMA spec- 
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FIG. 5: Low-frequency sector of the LMA impurity spectra, 
7rAoD(Cj) vs. u/Ao, for bandwidth D = 10~*Ao, i.e., u>o — 
O.OO8A0 = SOD, and interaction strengths U = 0.2Ao (top), 
U — 4Ao (middle), and U = 10 Aq (bottom). Discrete levels 
are represented by vertical lines. The Hubbard satellites at 
\io\ ~ f//2 lie outside the plot range. 



tra for the same parameters, along with an additional 
third calculation for U = IOAq — which is in the true 
strong coupling regime difficult to reach for the Lanczos 
method. 

Hofstetter and Kehrein's Lanczos spectrum for U — 
O.2A0 is dominated by Hubbard satellites at energies 
as high as |w| ~ ^ = 0.1 Aq = lOOOD which are ex- 
cluded from their plots.'^Slfn the physically more relevant 
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low-energy sector, their spectrum may be divided in two 
parts: (i) a central continuum for \uj\ < D = 10~''Ao 
which in shape is very similar to the Fermi-liquid band 
of the noninteracting limit; and (ii) two sets of narrow 
peaks on the antiferromagnetic scale O ( J ~ 12. 7D), the 
first located at \lu\ ~ 0.0018 - 0.002lAo (i.e., 18 - 21D) 
and the second at |w| ~ 0.0025-0.0026An (i.e., 25-26D). 
Although it is not clear how these sets will evolve with 
a larger number of sites in the Lanczos calculations, a 
structure with at least two features will most likely pre- 
vail in this region. 

The corresponding LMA results, displayed in the top 
panel of Fig. [5] reproduce — at least qualitatively — all 
characteristics of the Lanczos spectra: excellent agree- 
ment is found on the highest energy scale (off plot-range 
in the graph), governed by the Hubbard satellites at 
Lu ~ ±0.1 Aq with almost 97% of the total spectral in- 
tensity; good agreement is also observed on the low- 
est energy scale D, occupied by the central Fermi-liquid 
continuum, of net weight ©(ZJ/Aq) ~ 0.01%, and es- 
sentially unrenormalized from the noninteracting- limit, 
since, throughout the whole continuur n, th e interaction 
self-energies Tia{oj) ~ O (C/^/Aq) [eq. (39)] are weak in 
comparison to the hybridization A(a;) ^ ©(Aq), and 
have thus very minor influence on the quasi-particle prop- 
erties. Nevertheless, conceptually much simpler 2PT and 
UHF are similarly successful on the two latter energy 
scales. 

On the intermediate energy scale, J ~ 12.7-D, by con- 
trast, where UHF and 2PT both predict a single set 
of levels (see Sec. IV B), solely the LMA produces a 



rich structure which qualitatively resembles the Lanczos 
spectra.!^ The dominant features are two sets of poles at 
w = ±wi^2 [eq. (|4l|]: the first, with 2qi = 2.4% of the to- 
tal spectral intensity and located at = |a;i| ~ 1.28 J ~ 
O.OOI6A0 = 16Z?, corresponds to shifted "orbital rem- 
nants"; the second, with 2q2 — 0.6% of the net spectral 
weight and situated at \uj\ = \uj2\ — 0.78 J ~ O.OOIAq = 
lOD, is the novel feature related to resonant spin-flip 
processes (which has been discussed in the previous 
paragraph). Additionally, of the various sets of single- 
particle bands arising from the self-energy continua (see 
Sec. IIIB[ ) , only the two most prominent are visible in the 
graph: the first, at > ~ | J ~ 0.0019Ao = 19D 
and of width Z?, stems from spin-flipping impurity scat- 
tering of the metallic host electrons close to the Fermi 
level — a process driving the Kondo physics in the strong- 
coupling limit U ^ 4Ao (see Sec. IV D); a second, almost 



imperceptible set of continua of width 2D is located at 
> w< ~ J/2 ~ 0.000 64Ao = 6 AD and involves cou- 
pling the UHF orbitals to the low-lying spin-flip continua. 

In summary, contrary to UHF and 2PT, the LMA re- 
produces the main aspects of the Lanczos spectra also 
on the intermediate energy scale, but somewhat under- 
estimates the antiferromagnetic exchange J. A better 
quantitative match would require a renormalized spin 
exchange constant J' ~ 1.6 J, or equivalently a Coulomb 
coupling U' screened by the same factor for the spin chan- 



nel only. 

For U — 4Ao, not yet in the strong coupling regime, 
both the LMA and the Lanczos-determined spectra con- 
sist only of two appreciable contributions: the Hubbard 
satellites at w ~ ±2Ao on the high energy scale, and, on 
the other extreme, the low-energy Fermi-liquid contin- 
uum of width O {D) surrounding the Fermi level, a; = 0. 
The latter contribution, with its emergent central Kondo 
resonance, is plotted for the LMA in the middle panel 
of Fig. [5] The corresponding Lanczos graph by Hof- 
stetter and Kehrein confirms this scenario, albeit with 
a slightly less developed Kondo resonance, and sharper 
peaks at the band edges uj = ±1?. The overall agree- 
ment of the spectra is good, but could again be improved 
by renormalizing the LMA spin exchange, even though 
the mismatch might also partly stem from the artificial 
Lorcntzian broadening or the small number of sites used 
in the Lanczos calculations. 

Finally, the Fermi-liquid continuum for U — IOAq, 
plotted in the last panel of Fig. [5] illustrates the exponen- 
tial narrowing of the Kondo resonance with interaction 
strength which is one of the hallmarks of the strong cou- 
pling regime to be discussed in the next section. 



D. Strong coupling: U ^ 4Ao 

In an AIM with a narrow metallic host, the single- 
particle spectra suggest the following phenomenological 
definition of the strong coupling regime: an interaction 
strength belongs to the latter if (i) the metallic band at 
the Fermi level is dominated by a narrow central res- 
onance whose shape is almost identical to Kondo reso- 
nances belonging to even larger interaction strengths; (ii) 
the Hubbard satellites appear as sharply peaked continua 
instead of discrete levels; and (iii) the spectrum contains 
no other visible features. 

For U = 6A0 and 12 Aq, the LMA single-particle spec- 
tra do indeed comply with the last two of the above 
conditions. Whether however the shape of the Kondo 
resonance in the low-frequency sector of the spectra — 
displayed in the bottom row of Fig. [6] — is already scaling 
invariant is less obvious, but can be worked out in anal- 
ogy to the AIM with an infinitely wide metallic hostP^l^Il 
for U ^ 4Ao, and far from the band edges, <^ D, 
the contributions to Gct(w) = [a; — A(cj) — So.(i^)]^^ of 
both, uj and A-h^lo) ~ 2AoUj/nD, are negligible com- 
pared to the remaining two terms, given by Ai(cj) ~ Aq, 
and the interaction self-energy 'S^{uj). In the frequency 
range considered here, the latter is proportional to Aq 
and scales in terms of the single variable a) = lo/lo^, 
with uJm{U) ~ i:>exp[-7r[//8Ao] [see eq. (|40|]. This, in 
turn, leads to scaling invariance for the part of the single- 
particle continuum closest to the Fermi level, \uj\ D, 
i.e., for the central Kondo (or Abrikosov-Suhl) resonance. 



The limiting curve for U/Aq — >■ (X), plotted with red 
dashes in Fig. [7j is identical for any symmetric AIM 
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FIG. 6: Low-frequency impurity spectra, ■kAoD^lj) vs. u/Ao, 
for UHF (top), 2PT (middle), and LMA (bottom), for band- 
width D — O.OlAo and interaction strengths U = 6Ao (left 
column) and U = 12 Ao (right column). The high-energy 
Hubbard satellites are far off plot range. 



Fermi level as a result of an incipient transition towards 
an insulating atomic-limit state; 2PT, on the other hand, 
correctly captures the persistent metallic character of the 
system - manifest in the Friedel sum rule pinning of the 
spectra at the Fermi level - but misses the exponential 
narrovifing and the non-Lorcntzian shape of the Kondo 
resonance. 

Relatedly, on the high energy scale ^ U, the LMA 
correctly predicts sharply peaked Hubbard bands instead 
of genuine levels, ensuing from the absorption of the Hub- 
bard "levels" at |a;| ~ ^-|-^ by their former "sidebands" 
at f < |w| < f -f 2D; again, UHF and 2PT both fail to 
catch this many-body broadening effect. 

Finally, aU methods — UHF, 2PT, LMA, and also 
Lanczos — agree in the large spectral gap ranging from 
the low-energy continuum of width D to the high-energy 
Hubbard satellites. Especially the orbital remnants, 
which in weak and moderate coupling are found inside 
this gap, are now missing. A more rigorous analysis 
shows that they in fact still exist, occurring exponentially 
close to the band edges |a;| = D, but are undetectable be- 
cause their spectral weight vanishes exponentially. 
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FIG. 7: Scaling of the Kondo resonance in the LMA impurity 
spectra, nAoD^uj) as a function of u/ujm, with uj^ the (17- 
dependent) spin-flip or Kondo energy, for bandwidth D = 
O.OlAo and interaction strengths U/Aq — 4, 6, 8, 12 (solid 
lines from top to bottom), and in the limit U/Ao — > oo (red 
dashed line). Full discussion in text. 



with a metallic host. It comprises two parts, the first 
being the narrow central region between the cusps at 
\uj\ ~ Win, within which the line shape mostly carries the 
Lorentzian signature of Landau quasi-particles of weight 
Z = {1— [dReT,{uj) /duj]u,=o}~^ , while the second region, 
outside the cusps, is ch aract erized by long logarithmically 
decaying spectral tails.li^ESl (The cusps in the spectrum 
at |w| ~ Wm are artefacts stemming from the RPA-like 
structure of the polarization propagator and can be re- 
moved by a more realistic ansatz for the latter,'^ produc- 
ing LMA line shapes which then excellently agree with 
corresponding NRG data.l^ 

The graphs in the first two rows of Fig. [6] show clearly 
that such subtle effects are out of reach for the other 
two methods studied here: UHF, on the one hand, sup- 
presses for U/Aq — oo all spectral intensity close to the 



V. CONCLUSIONS 

In this article, a symmetric AIM with a narrow metal- 
lic host has been studied with three different theoretical 
approaches: UHF, 2PT, and the LMA. 

The persistent metallic character of the system - 
which renders diagrammatic perturbation expansions 
(like 2PT) viable in the first place - entails that its dy- 
namics hinges on the impurity Coulomb repulsion U and 
the noninteracting Green function g{io). The latter, in 
turn, depends solely on the hybridization A(w), which, 
for the present case of a narrow host band, is mainly de- 
termined by its width D and strength at the Fermi level 
Aq. Despite this reduced set of ultimately three param- 
eters - Aq, D, and U - the single-particle spectra are 
found to be rich, in particular in the physically relevant 
low-energy sector close to the Fermi surface, indicating 
the competition of various physical processes and their 
associated energy scales. These encompass Fermi-liquid 
behaviour within the low-energy continuum of width D, 
molecular orbital formation related to the bonding energy 
Wo ~ \/ AqD, antiferromagnetic phenomena driven by the 
exchange coupling J = 4ujq/U, and finally the Kondo ef- 
fect with its magnetic scale ~ Z?exp[— ttCZ/SAq]. 

Up to moderate interaction strengths, U < Aq, the 
low-energy physics is dominated by the integrated hop- 
ping between host and impurity, Uq — - J du!Ai{uj), and 
the antiferromagnetic exchange, J, which are both in- 
sensitive to the details of the hybridization function. The 
precise form of A(a;) therefore only determines the essen- 
tially unperturbed metallic band surrounding the Fermi 
level. In the opposite limit of large interactions, U ^ Aq, 
the low-energy physics is governed by the Kondo effect. 
This leads to a spectral scaling, in terms of w/wm, which. 
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again, does not depend on any details of the hybridiza- 
tion function. 

The LMA has been shown to produce meaningful re- 
sults over the entire range of interaction strengths. For 
small interactions, U ujq, it predicts a vanishing im- 
purity moment and connects smoothly to 2PT, yielding 
single-particle spectra carrying the signature of prevail- 
ing orbital physics. 

For moderate coupling strengths, ujq ^ U < Aq (or 
alternatively D <^ J <^ uiq), in addition to the high- 
energy Hubbard levels and the essentially unrenormal- 
ized Fermi-liquid continuum on the lowest energy scale, 
the LMA produces rich spectra on the J-scale: here, 
two pairs of poles along with several accompanying side- 
bands can be observed. Similar structures were found in 
corresponding Lanczos-determined spectra by Hofstetter 
and Kehrein.l^ Their main contributions can be rational- 
ized as follows: one pair of poles may be considered as 
remnants of the molecular bonding and anti-bonding or- 
bitals which similarly occur in UHF and 2PT; the other 
pair - which misses in the latter approaches and also in 
other state-of-the art techniques like, e.g., slave-boson 
based method^^ - is a somewhat unexpected feature 
and arises due to resonant collective spin-flip processes 
between the impurity and the UHF orbital remnants. 
Similar spectral contributions, in the form of sharp reso- 
nances at the inner band edges of the Hubbard satellites, 
have been observed numerically in the metallic phase 
close to the Mott transition, which occurs in the infinite- 
dimensional Hubbard modehf^Sisni These features have 
been dubbed antipolarons by Karski et al. who suspect 
them, mainly on energetic grounds, to take their ori- 
gin in bonding/antibonding phenomena between heavy 
quasiparticles and collective spin excitations.^ Although 
the same authors admit in a subsequent publicatiori^SI 
that this "complex composite excitation is not yet un- 
derstood," their original idea is supported by the above 
interpretation of the corresponding AIM features. Karski 
et al.'s results for the Hubbard model clearly indicate 
that sharp features will prevail in the presence of host or 



bath correlations, although their aspect may naturally be 
altered, especially if they coincide with another band. 

In the limit of large interaction strengths, U ^ Aq, 
of the three investigated methods only the LMA cap- 
tures simultaneously the many-body broadening of the 
high-energy Hubbard satellites and the relevant low- 
energy phenomena embodied in the exponentially nar- 
rowing Abrikosov-Suhl or Kondo res onanc e with its dis- 
tinct logarithmically decaying wings.!^^!^ Moreover, the 
observed Kondo resonance has been shown to possess the 
universal shape and scaling properties characteristic for 
symmetric ATMs with a metallic host. 

The fairly accurate analytic solutions obtained within 
the framework of the LMA suggest a classification into 
two regimes, which cover almost the entire range of in- 
teractions: the first englobes weak and moderate inter- 
actions, U < Aq, and its Fermi- liquid properties are pri- 
marily inherited from the noninteracting system, while 
its orbital contributions follow directly from the two-site 
approximation where the narrow host is treated as a sin- 
gle site; in marked contrast to the latter is the second 
regime, suitable for large interactions, U ^ Ag, where 
the narrow host band behaves as infinitely wide in com- 
parison to the exponentially small Kondo energy. 
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